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ABSTRACT 


Direct estimates of surface radiative fluxes that resolve regional and weather-scale variability 
over the whole globe with reasonable accuracy have only become possible with the advent of 
extensive global, mostly satellite, datasets within the past couple of decades. The accuracy of these 
fluxes, estimated to be about 10-15 W/m 2 , is largely limited by the accuracy of the input datasets. 
The leading uncertainties in the surface fluxes are no longer predominantly induced by clouds but 
are now as much associated with uncertainties in the surface and near-surface atmospheric 
properties. This study presents a fuller, more quantitative evaluation of the uncertainties for the 
surface albedo and emissivity and surface skin temperatures by comparingthe main available global 
datasets from the Moderate-Resolution Imaging Spectroradiometer product, the NASA Global 
Energy and Water Cycle Experiment Surface Radiation Budget project, the European Centre for 
Medium-Range Weather Forecasts, the National Aeronautics and Space Administration, the 
National Centers for Environmental Prediction, the International Satellite Cloud Climatology Project 
(ISCCP), the Laboratoire de Meteorologie Dynamique, NOAA/NASA Pathfinder Advanced Very 
High Resolution Radiometer project, NOAA Optimum Interpolation Sea Surface Temperature 
Analysis and the Tropical Rainfall Measuring Mission (TRMM) Microwave Image project. The 
datasets are, in practice, treated as an ensemble of realizations of the actual climate such that their 
differences represent an estimate of the uncertainty in their measurements because we do not possess 
global “truth” datasets for these quantities. The results are globally representative and may be taken 
as a generalization of our previous ISCCP-based uncertainty estimates for the input datasets. 

Surface properties have the primary role in detennining the surface upward shortwave (SW) 
and longwave (LW) flux. From this study, the following conclusions are obtained. Although land 


2 



surface albedos in the near-infrared remain poorly constrained (highly uncertain), they do not cause 
too much error in total surface SW fluxes; the more subtle regional and seasonal variations 
associated with vegetation and snow are still in doubt. The uncertainty of the broadband black-sky 
SW albedo for land surface from this study is about 7%, which can easily induce 5-10 W/nr 
uncertainty in (upwelling) surface SW flux estimates. Even though available surface (broadband) 
LW emissivity datasets differ significantly (3% - 5% uncertainty), this disagreement is confined to 
wavelengths > 20 pm so that there is little practical effect (1-3 W/m 2 ) on the surface upwelling LW 
fluxes. The surface skin temperature is one of two leading factors that cause problems with surface 
LW fluxes. Even though the differences among the various datasets are generally only 2-4 K, this 
can easily cause 10-15 W/nr uncertainty in calculated surface (upwelling) LW fluxes. 

Significant improvements could be obtained for surface LW flux calculations by improving 
the retrievals of (in order of decreasing importance): ( 1 ) surface skin temperature, (2) surface air and 
near-surface-layertemperature, (3) column precipitable wateramount and (4) broadband emissivity. 
And for surface SW fluxes, improvements could be obtained (excluding improved cloud treatment) 
by improving the retrievals of (1) aerosols (from our sensitivity studies but not discussed in this 
work), and (2) surface (black-sky) albedo, of which, NIR part of the spectrum has much larger 
uncertainty. 
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1. Introduction 


Analyses of global satellite products have resulted in the production of several extensive and 
detailed surface radiative flux datasets [Rossow and Lads, 1990; Darnell and Staylor, 1992; 
Whitlock et al., 1995; Rossow and Zhang, 1995; Stackhouse et al., 2004; Zhang et al., 2004]. 
However, the usefulness of these flux products for monitoring long-term variations of global 
radiation budgets and climate changes, as well as in other applications such as biological and 
oceanic modeling, is dependent on two separate issues: (A) the current uncertainty (or reliability) 
that can be estimated by evaluation of the calculated radiative fluxes against more direct 
measurements and (B) how much and by what means may the flux products be improved. To address 
(A), Zhang et al., [1995 and 2004] have compared their flux products against surface (and the top 
of atmosphere) flux observations to estimate their uncertainty. Their 2004 estimate is about 10-15 
W/m 2 uncertainty for surface fluxes. For (B), they have conducted all the important sensitivity 
studies by varying the input datasets (and the radiative transfer model parameters) to provide 
quantitative uncertainty estimates of the impacts of realistically-assumed uncertainties of the input 
datasets. The “realistically-assumed” uncertainties for the input datasets were based on earlier 
estimates, mainly from the International Satellite Cloud Climatology Project (ISCCP) D1 [Rossow 
and Schiffer 1991]. Some of their conclusions are summarized here. (1) Although clouds have, for 
a long time, been highlighted as the major source of uncertainty in both planetary and surface 
radiation budgets, the advent of extensive cloud datasets has reduced this source of uncertainty. (2) 
As a result, other uncertainty sources have now become relatively more important and comparable 
to the uncertainty from cloud datasets. (3) For surface radiative flux estimates in particular, the 
accuracy is now largely limited by two aspects of Earth observations that have long been assumed 
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to be adequately observed, especially during the satellite era, namely, the near-surface atmospheric 
radiative properties (temperature and humidity) that are usually obtained from atmospheric profiling 
instruments and the surface radiative properties (surface skin temperature, solar albedo and infrared 
emissivity). 

Zhang et al. [2006] have presented a fuller, more quantitative evaluation of the uncertainties 
for the near-surface air temperature and humidity by comparing the main available global datasets 
that are treated as an ensemble of realizations of the actual climate such that their differences 
represent an estimate of the uncertainty in their values. Here we follow the same methodology to 
make such comparisons for the surface properties (Section 2). In Section 3 we briefly summarize 
the uncertainties of the surface radiation budget induced by the uncertainties of these input datasets 
and discuss some implications for understanding andmodeling the surface-atmosphere interactions. 


2. Surface Radiative Properties 

The five major global datasets (three from reanalysis and two from satellite-retrieval), used 
in the comparison for near-surface atmospheric properties in our companion paper [Zhang et al. 
2006], also supply other parameters such as surface skin temperature and surface albedo, which are 
to be compared in this work (with some additional datasets). They are, (1) the first version of the 
ECMWF reanalysis for a 15-year period (1979-1993, called ERA15 or ER in our figures/tables) 
[Gibson et al., 1999], (2) the first version of the NASA Data Assimilation Office (DAO) Goddard 
Earth Observing System (GEOS) reanalysis fora 15-year period (March 1980 to Februaryl995, 
called GEOS-1 or GE in our figures/tables) [Takacs, et al., 1994], and (3) the second version of the 
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National Centers for Environmental Prediction reanalysis for a 51-year period (1948-1998, called 
NCEP or NC in our figures/tables) [Kistler et al., 2001], (4) a version of TOVS from the 
International Satellite Cloud Climatology Project (ISCCP, thereafter ISCCP-TOVS) [Rossow et al., 
1996] that is combined with the modified surface skin and air temperature from ISCCP-FD [Zhang 
et al., 2004], called ISCCP-FD-TOVS, shortened as FDTV (called TV in some figures/tables, and 
(5) 31 dataset from the Laboratoire de Meteorologie Dynamique (LMD) group, who has developed 
a new analysis that produced temperature and humidity profiles dataset from TOVS (Path B) 
[consisting of the High resolution Infrared Radiation Sounder (HIRS), the Microwave Sounding Unit 
(MSU) and the Stratospheric Sounding Unit (SSU)] based on Improved Initialization Inversion 
(therefore, 31) retrieval algorithm [see Scott et al,, 1999]. For more details of these five datasets, see 
Zhang et al., [2006]. To simplify the presentation, unless otherwise indicated, we restrict the 
comparisons to monthly means and, if available, monthly-hourly means for every sixth hour (UTC 
= 0\ 6 h , 12 h and 18 h ), for January, April, July and October of 1992. Except for some long-term 
variations in the character of these datasets that we note below, this year is taken to be typical of 
these products. 

2.1. Albedo 

The surface albedo, as defined by the ratio of the reflected to incident flux density at surface 
and reported from measurements or calculations, is usually the apparent albedo that is governed by 
the radiative interactions between atmospheric multiple scattering and absorption and surface 
reflection and absorption under specific atmospheric conditions (including clouds) and is also solar- 
zenith-angle dependent. As a result, it is highly variable in time because of the rapid variation of 
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atmospheric conditions and the slower variations of surface conditions. For theoretical and practical 
uses, particularly in radiation modeling that can account for multiple scattering, there has long been 
used the physical concept of a surface albedo without atmospheric (including clouds) effects. 
Though, this quantity is sometimes referred to “surface albedo” without clarification by many 
authors, a number of authors have actually distinguished it from the apparent albedo by using 
different names, e.g., “black-sky albedo” (BSA, as chosen in this study) [Jin et al., 2003b], “true 
albedo” [Henderson-Sellers and Wilson, 1983], and “inherent albedo” [Liang et al, 1999]. Thus, 
BSA is the intrinsic property of the surface material (but still varies with solar zenith angle), 
independent of the properties of any overlying layer (atmosphere). Confusion arises when people 
use the word “albedo” and neglect the difference between the apparent and black-sky albedo. 
Henderson-Sellers and Wilson [1983] actually pointed out: “Atmospheric scatter tends to increase 
observed clear-sky system albedo compared with true surface albedo.” Sometimes, people just 
simply and approximately take “clear-sky albedo” (which is a special case of the apparent albedo, 
but still is very sensitive to the aerosol loading and humidity profile in the atmosphere) as BSA. 
Note BSA is the apparent albedo for the moon’s surface or for other no-atmosphere planets. To 
derive surface BSA from actual measurements on the Earth, the measured (apparent) albedo values 
must be “corrected” to eliminate all atmospheric multiple scattering and absorbing effects, including 
shifts of the spectrum of illumination and the contribution from multiply-scattered radiation. BSA 
is widely used because it is BSA that is specified in a radiative transfer model for calculating 
shortwave (SW) fluxes. In the following text, we use “albedo” alone when it is for general purposes, 
meaning that it can be any one of the different albedos or when the precise meaning is not known, 
while we use BSA or clear-sky albedo in a precise way. Similarly, the longwave (LW, terrestrial 
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radiation wavelengths) surface emissivity can be determined with or without atmospheric effects 
(no-atmosphere or “black-sky emissivity”, section 2.2). 

A. Datasets 

Since the advent of satellites, various global, broadband, land surface albedo products have 
been derived either from satellite measurements alone [e.g., Staylor and Wilber, 1990, and Li and 
Garand, 1994, both from the Earth Radiation Budget Experiment (ERBE); Jin, et al., 2003b, from 
the Moderate -Resolution Imaging Spectroradiometer (MODIS)], or from a combination of the 
surface spectral reflectivity properties for different surface types, compiled from surface, airborne, 
satellite and laboratory measurements, and classification of each location by vegetation-soil types 
[e.g., Csiszar and Gutman, 1999; Zhang et al., 1995]. The latter approach is used in virtually all 
weather forecast and climate models. Both methods (or some combination) are used by various 
authors in their calculations of surface and/or the top of the atmosphere (TOA) SW fluxes [e.g., 
Rossow andLacis, 1 990; Darnell and Staylor, 1992; Pinker and Laszlo, 1992; Whitlock, etal., 1995; 
Zhang et al., 1995; Chen and Roeckner, 1996; Gupta, et al., 2001 ; Zhang et al., 2004]. 

ERA 15 used an albedo map [ Gibson , et al., 1999] determined from a background yearly 
climatology with fixed values of 0.55 for sea ice and 0.07 for open water. For land the albedo varies 
between 0.07 and 0.80 whenever snow occurs but the albedo for snow-free vegetation is constant 
all year long, based on surface vegetation classification. GEOS- 1 specifies its surface albedos from 
Posey and Clapp [1964] with some modifications for land and fixed values of 0.07 and 0.80 for open 
water and sea ice, respectively. For snow cover, albedo is usually greater than 0.4 [ Takacs , et al., 
1994]. The NCEP surface albedo comes from a climatology [ Kanamitsu , 1989; Kalnay, 1996]: the 
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current land albedo is from Matthews [1985] with modification if there is snow cover (0.75 is 
assumed poleward of 70° latitude for permanent snow and 0.60 is assumed equatorward of 70° if 
snow depth is at least 0.01 m), sea ice albedo ranges from 0.45 - 0.65 and 0.65 - 0.8 for snow-free 
and snow-covered, respectively, and open water albedo is solar-zenith-angle dependent but no 
specifics could be found. NCEP albedos do not depend on the spectral interval (see 
ftp://ftpprd.ncep.noaa.gov/pub/cpc/wd51we/reanal/random_notes/model). Since all of the above 
albedo models are relatively simple, we do not use them for this work. 

The NASA Global Energy and Water Cycle Experiment (GEWEX) Surface Radiation 
Budget (SRB) Project has developed a 12-year/ 148-month global dataset for surface SW and LW 
fluxes on 1° X 1° grid, called SRB release 2 [ Stackhouse et ah, 2001, 2004] (now extended to 22 
yrs). The SW algorithm used for SRB is from Pinker and Laszlo [1992] in which land and ocean 
albedo values, as a function of solar zenith angle and wavelength, are clear-sky surface albedos 
based on Briegleb et al. [1986]. We call this version of albedo SRB-ALB. 

The NASA Goddard Institute for Space Studies (GISS) GCM Model II [Hansen, etal., 1983] 
uses a global, seasonally varying land surface albedo dataset that was constructed by combining a 
land surface type/vegetation classification from the work of Matthews [see Matthews, 1985] with 
field and laboratory albedo measurements for different rock, soil and vegetation types. Snow is 
treated as an added component that modifies the surface albedo accordingly. Both snow and ice 
albedos are also (implicitly) solar-zenith-angle dependent. The spectrally and solar-zenith-angle 
dependent ocean albedo is based on calculation of Fresnel reflection from a surface wave-slope 
distribution as a function of wind speed [ Cox and Munk, 1956], The 1983-version of the albedo 
values are for two wavelength ranges: visible (VIS) from 0.2 to 0.7 qm and near-infrared (NIR) from 
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0.7 to 5.0 |j,m, together coveringthe whole solar spectrum The revised 200 1 (cunent) version of this 
product, which we will call GCM-ALB, has albedo values for six wavelength ranges: VIS and five 
covering the NIR. The revised ocean albedo includes the effects of foam and hydrosols. For all of 
the work using the GISS model, the default wind speed= 2 m/s [for a more detailed description, see 
Zhang et al., 2004] . The version used here for albedo values has a spatial resolution of 280 km and 
is “climatological” in that it represents no particular year (except for snow/ice areas). 

Based on the 1983 version of the GISS GCM, Zhang et al. [1995] produced a broadband 
albedo product by combining the surface visible albedo taken directly from visible reflectances (at 
wavelength about 0.6 |am) in the ISCCP cloud products [Rossow and Garden, 1993] and the GCM’s 
ratio of the NIR to visible albedo. The land-type-dependent visible albedo is replaced by a 3 -year 
seasonal average of the ISCCP-C1 visible reflectances and the NIR albedos adjusted by regression 
with ERBE’s TOA clear-sky SW fluxes to minimize systematic errors. This product has a spatial 
resolution of 280 km and represents albedo values for particular time periods (1985 - 1988). We 
refer to this dataset as ISCCP2-ALB (where the “2" refers to the number of spectral bands). 

Similarly, we have produced a updated version of the previous product by employing the 5 
spectral ratios of the 5 NIR wavelengths to visible albedo from the 2001 GISS GCM and the new 
ISCCP-D1 surface visible reflectances (now corrected for aerosols using a climatology). We refer 
to this dataset as ISCCP6-ALB, which is used in our new radiative transfer flux product ISCCP-FD 
[Zhang et al., 2004]. 

The MODIS team (of Boston University) has produced a land surface albedo dataset at 
several spatial resolutions (from 5 km) for 16-day periods since February 2000 based on a three- 
parameter semi-empirical RossThick-LiSparse-Reciprocal (RTLSR) bidirectional reflectance 
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distribution function (BRDF) model to characterize the anisotropic reflectivity of the land surface 
[Jin et al., 2003a and b]. In their albedo retrieval irradiance model (MODTRAN), twenty surface 
reflectance spectra, four visibility values for different aerosol loadings, five default aerosol models, 
twelve different water vapor profiles from the Tropical Ocean Global Atmosphere Coupled Ocean- 
Atmosphere Response Experiment project but with adjustments, and five types of default clouds 
(essentially some climatological atmosphere) are used [Liang et al., 1999]. The snow/ice flag is 
turned on if the majority of observations are indicated as snow/ice during a 16 day period, and only 
those observations are used for BRDF and albedo snow inversion [Schaaf et al., 2002] . The surface 
cover types are assumed Lambertian. In this study, we use a dataset provided by Crystal B. Schaaf 
at http://duckwater.bu.edu/brdf_albedo/albedol6/params.htm from the original version 4 product 
(MOD43C2 V004 16-day L3 Global 0.05° CMG, as of 2004). The original BRDF parameters for 
0.05° X 0.05° cells are used to calculate BSA for three wavelength ranges: visible (0.3-0. 7 |im), NIR 
(0. 7-5.0 |j.m) and total SW (0. 3-5.0 |j.m) using a polynomial function of solar zenith angle (supplied 
byMODISteam, seeLucht et al., [2000])witha 16-day-mean solarzenith angle (representing a 16- 
day period “mean” BSA). The fine-scale albedo values are then (linearly) averaged to our standard 
280-km equal-area-map grid cells when > 90% of the 0.05° cells have values within a box. This 
albedo dataset, called MODIS-ALB, is for land only. Using the fonnulae provided by theMODIS 
team, we have produced MODIS BSA for four 16-day periods (beginning dates are 6 April, 1 1 July 
and 29 September 2000, and 17 January 2001) for comparison with matched values from ISCCP6- 
ALB for visible, NIR and total SW BSA, respectively. In addition, the MODIS BRDF parameters 
for one 1 6-day period (11-24 July 2000, representing all of July) and two 16-day periods ( 1 January 
to 1 February 200 1 ) are used to calculate 3-hourly surface visible and NIR BSA values for these two 
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months that are then used in our current radiation model [Zhang et al., 2004] in place of ISCCP6- 
ALB to calculate ISCCP-FD-like MODIS (thereafter, MODIS-FD) 3-hourly TOA/surface clear-sky 
albedoes/fluxes for July 2000 and January 2001. The monthly means are compared with matching 
values from ISCCP-FD and the Clouds and the Earth’s Radiant Energy System (CERES) ERB E-like 
datasets [ Wielicki et al., 1996; see also Zhang et al., 2004]. 

There are also various regional albedo products, e.g., the Meterosat Surface Albedo based 
on two adjacent geostationary satellites Meteosat-7 and 5 (currently operated by EUMETSAT). The 
production relies on a daily accumulation of the satellite observations acquired under more than 
fifteen different illumination conditions to assess the scattering properties of the surface and the 
atmosphere with some assumptions, and estimates the directional hemispherical reflectance values, 
for visible band (~ 0.5-0. 9) that is corrected for atmospheric effects (i.e., BSA), which may be 
linearly related to broadband albedo. The product is planned for retrievals of albedo values for the 
whole Meteosat record (1982 to current) on a daily basis [Govaerts et al., 2004]. 

B. Comparison of Broadband Albedos 

Table la shows the global mean (rms = rms with bias removed, this will be used for all rms 
values throughout the paper) of regional differences for monthly-mean clear-sky surface albedo 
between ISCCP6-ALB and SRB-ALB for January, April, July and October, 1992 for four broad 
surface types: snow-free land, ice-free water, snow-covered land and sea ice, where snow/ice cells 
are defined by coverage >95% of a 280-km equal-area box on all the days of a month based on the 
ISCCP snow/ice data. The lastrow of the table summarizes the average results for the four seasonal 
months. In general, the two albedo sets are in good agreement for snow-free land and ice-free water: 
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the mean (rms) differences (ISCCP6 minus SRB-ALB) are -2.0% (4.6%) and -0.3% (3.0%), 
respectively. But for snow-covered land and sea ice, the differences become 11.2% (6.5%) and 
11.0% (8.3%), respectively. We have also compared these products for seven representative 
vegetation types: (1) Amazon rainforest, (2) Sahara Desert, (3) snow-free grassland, (4) snow-free 
deciduous forests, (5) snow-covered grassland, (6) land ice (with or without snow) and (7) sea ice. 
These regions are defined in the same way as above except that the snow-free vegetated regions are 
defined by > 66 % coverage by the type of vegetation for a 280-km equal-area grid cell. Table lb 
shows the comparison between ISCCP6-ALB and SRB-ALB for the four monthly mean clear-sky 
albedos for the seven surface types. The largest biases (magnitude > 5%) appear in Amazon 
rainforests and snow-covered grassland: albedo values are more difficult to be detennined in these 
regions because of contamination by persistent clouds and snow, respectively. 

Table 1 c compares surface BSA for visible, N1R and total broadband SW between 1SCCP6- 
ALB and MOD1S-ALB for the four (seasonal) 16-day periods for snow-fitee land (defined by both 
snow maps being snow-free), snow-covered land (both maps have grid cells with > 75% snow) and 
all land (all available cells for both). For snow-free land, there is a good agreement in the visible 
band: mean (rms) regional differences (FD minus MOD1S) are 0.9% (3.5%). Disagreement in the 
N1R is much greater: mean (rms) regional differences are -19% (6.7%). The MOD1S retrieval 
algorithm has larger errors in theNIR than in the visible (. Jin et al., 2003a) because the radiances 
are smaller and the atmospheric absorption is much stronger. On the other hand, the N1R albedo 
values in 1SCCP6-ALB rely on a surface-classification-based GCM model, which is a long-tenn 
climatological surface representation based on Matthews [1985], but adjusted by regression with 
ERBE to reduce systematic errors. There is not sufficient evidence to decide which of these N1R 
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results is better than the other (see next paragraph). For snow-covered land, both the agreement and 
disagreement are smeared out and their mean (rms) regional differences become -3% (14%) and - 
6.8% (10%) for visible and NIR, respectively. 

The visual comparison for the seven representative vegetation types is shown in Figure 1 . 
Figures la to lc are scatter plots between ISCCP6 (=X) and MODIS (=Y) BSA for Amazon 
Rainforests, snow-free grassland and snow-free Deciduous Forests for visible, NIR and total SW, 
respectively, for the same four 16-day periods (on 280-km equal-area map), and Figures Id to If are 
their counterparts for snow-covered grassland, land ice, sea ice and Sahara Deserts. Because of 
sparse sampling in MODIS, there are only 7 and 6 grid cells with available albedo values from all 
the four periods for land ice and sea ice, and as a result, their comparisons are statistically 
insignificant. For the other five regions, generally speaking, the figures are more or less different 
from the mean comparison results shown in Table lc. The Sahara desert and snow-covered 
grasslands have the highest correlations (>0.71) for all the three bands, but large mean differences 
always appear in NIR (-22% and -10% in ISCCP6 minus MODIS, respectively). The least 
agreement appears in Amazon Forests (correlation coefficients are 0.46 and 0.41 for visible and 
NIR, respectively while their mean differences are 0.7% and -22%, respectively). Also, ISCCP6 and 
MODIS have better agreement in visible band than in NIR band as also shown in Table lc 

To obtain a global picture, figure 2a to 2c show the global maps of the BSA differences 
(ISCCP6 minus MODIS) for visible, NIR and total SW, respectively, for the 16-day period of July 
1 1 (beginning date). Again, the visible map shows that, for the majority of the world, the two dataset 
differ only within ± 5% (mean global difference = -0.5%). The largest discrepancy (up to -35% 
difference) appears in central Greenland (where MODIS has less data), possibly because of the 
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difficulty of identifying clear conditions over such bright regions. There is a narrow belt from 
England to the Northeast China for essentially intensively-cultivated lands [Matthews, 1982] where 
ISCCP6 is 5- 1 0% higher than MODIS . The NIR map shows large difference (mean global difference 
= -20%) and systematical regional-pattern dependence, indicating the two datasets have some 
fundamental difference in deriving the NIR albedos. The most prominent region is the North African 
deserts, where ISCCP6 is less than MODIS by > 30%. For the total broadband SW, the majority of 
land areas show a difference around -10% (also the global mean difference). 

Table Id shows the comparison of calculated monthly-mean surface clear-sky (apparent) 
albedo and upward SW between ISCCP-FD and MODIS-FD for snow-free, snow-covered and total 
land. The agreement between ISCCP-FD and MODIS-FD is very good showing that the large NIR 
BSA differences have far less effect on the total clear-sky (apparent) albedo and upwelling SW 
fluxes because the atmosphere absorbs most of the SW radiation in the NIR band. For all land, the 
mean (rms) regional differences between ISCCP and MODIS are only 0.14% (7%) and -2.5 (11) 
Wm 2 for the albedo and upwelling SW, respectively. 

Tables 2a and 2b show the comparisons of ERBE TOA clear-sky albedo with ISCCP-FD and 
SRB, respectively, for the four seasonal monthly means of 1986 when ERBE data is available. For 
snow- and ice-free lands, both TOA clear-sky albedos are in good agreement with ERBE: ISCCP 
is within 1% bias while SRB is within 2 %. This is not the case for snow-covered grassland, land 
and sea ice (land types 5 through 7), where ISCCP-FD exhibits a bias up to about 8% while the SRB 
bias is up to 13%. For snow-covered grassland (type 5), ISCCP-FD is biased high while SRB is 
biased low; but for land and sea ice, both ISCCP-FD and SRB are biased low. The reasons may lie 
partly in their fundamentally different analysis approaches and partly in the cloud contamination for 
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ERBE. One known problem for the high latitude results in SRB, which uses the ERBE angle models 
to determine fluxes from radiances [Pinker and Laszlo, 1992], is that these angle models are biased 
for large solar zenith angles [Suttles et al., 1992]. 

Although ISCCP-FD and MODIS have mean regional differences of -7.6% for surfaceBSA 
(based on four 16-day periods and all land areas), these differences produce TO A clear-sky albedo 
mean (rms) differences of only -0.31% (3.4%) (not shown) because TOA albedo is insensitive to 
surface NIR albedo, the visible band contributing most to the SW fluxes at TOA. Table 2c and 2d 
compare monthly mean TOA clear-sky albedos fromCERES (ERBE-like) with those from ISCCP- 
FD and MODIS-FD, respectively, for July 2000 and January 2001. The bias of ISCCP-FD with 
respect to CERES (-0.29%) is slightly smaller than for MODIS-FD (0.70%) for snow-free land, but 
the rms difference is slightly larger (2.4%) than for MODIS ( 1 .7%). For snow-covered land, ISCCP 
has a larger bias (4.6%) and rms difference (1 1%) than MODIS (3.0%, 9.8%). For all land, ISCCP- 
FD and MODIS have biases (nns) = 0.90% (5.7%) and 1.2% (4.5%), respectively. Both albedo 
products produce good agreement with the CERES albedos on average but ISCCP exhibits more 
scatter relative to CERES than MODIS does. These statistics do not provide strong evidence for 
choosing between these two albedo products. Figure 3 shows the scatter plots of TOA clear-sky 
upward SW fluxes from CERES against the results from ISCCP-FD and MODIS-FD for snow-free 
and snow-covered land areas. From figures 3a and 3b (snow-free land), we see that MODIS appears 
to perform slightly better than ISCCP: the points off the diagonal (> 1 10 Wm 2 ) in Figure 3a are 
located in Northern Africa and the Arabian Peninsula, where ISCCP is biased low. The fact that the 
pattern of disagreement for snow-covered land shown in figure 3c and 3d is the same for both 
ISCCP-FD and MODIS-FD suggests that this disagreement might be partly caused by the ERBE 
angle models for this surface type. 
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Table 3a (for snow-free landand ice-free ocean) and 3b (for snow-covered land and sea ice) 
compare TOA albedo and upwelling SW fluxes for clear-sky scenes with ERBE for ISCCP2-ALB, 
GCM-ALB, ISCCP6-ALB and SRB-ALB. The albedo and fluxes are calculated from their 
corresponding surface broadband BSA or clear-sky albedo, whichever applies. From the first three 
entries, we can see that adjusting the GCM albedo using the ISCCP visible reflectances for land 
improved the albedo/flux comparison substantially. Going to six wavelength ranges also provides 
some overall improvement over the original two- wavelength-range treatment. For snow- free land 
and ice-free water regions, the SRB-ALB results are comparable with ISCCP6-ALB, i.e., both have 
a good agreement with ERBE, but for snow-covered land and sea ice regions, SRB-ALB has a larger 
bias than ISCCP6-ALB. There are three notable conclusions. First, except for GCM-ALB, all of the 
other products are adjusted to agree with ERBE in some way, yet the remaining regional scatter in 
surface albedo for snow-free land or ice-free water is still several percent, and even larger for 
snow/ice surfaces. Second, although the reconstruction of surface albedo from a surface 
vegetation/land use classification produces a fairly accurate result as shown for NASA GISS GCM 
model (which can be compared with other major GCM results), use of the ISCCP visible reflectance 
data improves the result, which implies that there is significant regional (and temporal) variation in 
the albedo of surfaces with the same classification [cf. Matthews and Rossow, 1987]. Third, it is still 
possible to improve surface albedo maps by producing and using better satellite-based surface 
albedo datasets. We note that MODIS is the first imaging instrument to provide measurements at 
wavelengths covering the whole solar spectrum. 

Figure 4 shows the zonal mean TOA clear-sky albedo from ISCCP6-ALB (from ISCCP-FD), 
SRB-ALB and ERBE for July 1986 and January 1987. In most zones, the former two are 
comparable and exhibit biases with ERBE < 5% or so. In the winter hemispheres, at about 40°S- 
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55°S for July and 30°N-45°N for January, SRB performs better thanISCCP-FD, while in summer 
hemispheres, ISCCP usually performs better than SRB. In the polar regions, ISCCP-FD values are 
closer (especially in January) to ERBE than SRB. 

2.2. Infrared Emissivity 
A. Datasets 

Infrared emission from the surface depends on both temperature, which varies rapidly in time 
especially for land, and emissivity that varies with wavelength and surface properties in a complex 
way. Emissivity also plays an important role in determining surface temperature by remote sensing 
[e.g., Becker, 1987]. Because the thermal (LW) emissivity at the Earth’s surface is close to unity, 
most applications assume it to be unity [ Peixoto and Oort, 1992], which is also the case in Zhang 
et al. [1995]. As a result, there were no high-quality global surface emissivity datasets determined 
from satellite observations (31 produces emissivity at 50 Ghz, see also Prigent et al., [1997 and 
2000], but not for thermal radiation spectral range) or otherwise at the start of this work. The only 
global databases for this surface property (in terms of LW flux calculation) are constructed by 
combining field and laboratory measurements of the emissivities of various rock, soil and vegetation 
types, which are generally only measured for wavelengths < 20 qm, with a surface classifications 
and theoretical calculations (especially for water as in the GISS climate model, where the emissivity 
= 1 - albedo and the albedo is calculated using Fresnel reflection as explained before). 

Since July 2001, there has appeared several versions (still kept updating) of the MODIS- 
based high-spatial-resolution Land Surface T emperature product that also includes emissivity, which 
is for 6 wavelength bands, ranging from 3.66 to 12.17 qm [see e.g., Wan et al., 2004]. Because this 
product is only for several narrow bands and still being validated, we do not used it in this study. 
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Two versions datasets are examined here, one from the GISS GCM and used in ISCCP-FD 
with a spatial resolution of 280 km and emissivity values in 33 k-intervals covering the whole LW 
range (5 -200 pm) [cf. Lacis and Oinas, 1991] and one adopted from Wilber etal. [1999], used for 
the determination of surface radiative fluxes by the SRB and CERES with the Fu-Liou radiative 
transfer code [Fu andLiou, 1992], In deriving the broadband emissivity (5.0-100 pm), Wilber etal. 
[1999] used laboratory measurements for nine materials (including water) over a wavelength range 
of 4-16 pm and, by averaging 12 Fu-Liou band values (covering4.5 -°° pm)weighted bythe Planck 
function energy distribution, extrapolated them to broadband emisivities for 18 surface types 
(including water), representing all the possible surface types. The global broadband emissvity (for 
nadir viewing conditions) is then constructed based on the 1 8 types on a 1 0’ X 1 0’ global grid, which 
we regridded to 280-km equal-area map for our comparison. We refer to these two datasets as GCM- 
EMIS and SRB-EMIS, respectively. In the processing, Wilber et al. [1999] have also derived the 
12 Fu-Liou band-average emissivities as well as CERES window (8-12 pm) emissivities for the 18 
surface types. 

The ISCCP dataset retrieves land surface brightness temperatures (with emissivity assumed 
to be unity) from infrared radiances measured at ~ 1 1 pm wavelength in cloud-free scenes and 
corrected for atmospheric emission/absorption [cf. Rossow and Gander, 1993]. If accurate values 
of the land surface infrared emissivities at 1 1 pm are available, these data can be corrected to 
estimate surface skin temperatures (T s ). Hence, we also consider the available information about the 
“narrowband (11 pm)” or window region infrared emissivity. The version considered here also 
conies from the GISS GCM database (11. 1-11. 3 pm) and is compared with the CERES emissivity 
(8-12 pm). 
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B. Comparison of Broadband Emissivities 

Figure 5 shows a map of the differences of broadband emissivities between GCM-EMIS and 
SRB-EMIS. Values over land generally agree to within ±2%, but the coherent features exhibited in 
the figure suggest regional differences in classification. The GCM-EMIS values tend to be 
somewhat larger than the SRB-EMIS values for denser vegetation types (e.g., rainforests) and for 
some heavily cultivated locations (e.g., Europe and China, but not in Southern Africa), and 
somewhat smaller than SRB-EMIS for intennediate density grasslands (e.g, the Russian Steppes, 
central US) and most of middle/southern Africa (from dense forest to grass savanna), for 
mountainous terrains and even smaller in frozen areas (e.g., Greenland and Antarctic). 

Surprisingly, the largest differences are for open water (and ice/snow), ranging up to 7%. The 
SRB water emissivity neglects not only the strong spectral dependence [since its broadband 
emissivity is extrapolated from 1 2-band-averages over 4- 1 6 pm (weighted with the Planck function 
energy distribution)], but also zenith angle variations and windspeed effects. Averaging over all 
zenith angles produces an emissivity, even at “window” infrared wavelengths that is lower than the 
nadir value (usually measured in a laboratory). Wind roughening of the surface is equivalent to re- 
distributing the zenith angles, but this has only a weak effect on the angle-averaged emissivity. For 
a calm ocean surface, the emissivity (actually the reflectivity) can be accurately calculated using the 
Fresnel reflection formula [e.g., Lion, 2002], and under windy conditions, such calculations can be 
done with sufficiently high precision using a wave-slope distribution model as used in the GISS 
model. Figure 6a illustrates the variations of the GISS GCM emissivity values for water at 
wavelengths from 0.2 to 200 pm and surface wind speed = 3 m/s for cosine zenith angles 0.1, 0.5 
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and 1.0, respectively, and Figure 6b shows the emissivities at cosine zenith angle = 0.5 for wind 
speed = 3 and 7 m/s, respectively from Jacek Chowdhary (personal communication). From Fig. 6a, 
we see that the spectral difference in emissivity can be > 15% and the extrapolation of SRB’s 
emissivityto wavelengths beyond 16 pm tends to overestimate the broadband emissivity [5-100 pm 
as obtained by Wilber et al., 1999]. 

The difference of snow/ice surface emissivity comes from the use of different laboratory 
measurements for the GCM-EMIS and SRB-EMIS, the former adopted laboratory results from a 
variety of sources while the latter is basically from Salisbury and D ' Aria [1992], 

C. Comparison of Spectral Emissivities 

To convert the ISCCP surface (skin) brightness temperatures into physical temperatures, we 
use the narrowband emissivity at about 1 1 pm wavelength [see discussion about the trade off 
between correcting temperature and correcting emissivities in the work of Zhang et al, 2004]. 
Figure 7 shows the global difference map between the GCM-EMIS and SRB-EMIS narrowband 
values. Here the ocean values are in very good agreement while land and snow/ice locations 
generally agree within ±3%. The largest disagreements occur over arid and mountainous locations, 
including tundra: in all of these regions the GCM-EMIS values are significantly larger (> 1 0%) than 
SRB-EMIS. Examining the reported values in detail suggests that the latter dataset employs values 
consistent with bare rock and soil in theseregions, whereas the former assumes some effect of sparse 
vegetation for most of these regions. The narrowband emissivity displays similar dependence on 
zenith angle as shown in Figures 6a but the surface windspeed dependence is much weaker. To 
correct the ISCCP surface temperature retrievals, which are made predominantly at small zenith 
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angles (the average value is less than about 30°), we use a single value (0.985), a mean windspeed 
of 2 ms 1 and a simple approximation of the latitudinal dependence of reflected downwelling 
radiation for ISCCP-FD calculations [Zhang et al., 2004], In our sensitivity tests [see also Zhang et 
al., 2004], if wind speed changes from 2 to 10 ms" 1 , the emissivity change would induce a upwelling 
LW flux mean (rms) change by 0.5 (0.3) W/m 2 for global ocean, while a ±2% change (i.e., 4% of 
change) of emissivity would induce 2.1 (1.3) W/m 2 change for the whole globe. 

2.3 Comparison of Surface Skin Temperatures 

All the major datasets have an associated surface skin temperature (T s ) as shown in Zhang 
et al. [2006: Table 1], but the values from FDTV and 31 are retrieved from satellite radiances using 
their respective temperature/humidity profiles (the original FDTV, i.e., ISCCP-TOVS, T s values 
come from the ISCCP analysis that combines infrared radiances from imaging radiometers with the 
TOVS atmospheric profiles, whereas 31 uses the “window” channel radiances that are part of the 
temperature/humidity profile retrieval) and the other three are determined by the assimilation model 
as explained in Zhang et al. [2006: section 2.1]. 

Comparison of the monthly-mean T s maps from the five datasets (ERA 15, GEOS-l,NCEP, 
31 and FDTV) for the four seasonal months shows strongly similar pictures with spatial correlation 
coefficients for all possible pairs ^0.95 (0.980 overall). Ocean areas are better correlated (0.986) 
than land areas (0.975). The largest regional differences, > 40 and 35 K for land and ocean, 
respectively, between FDTV and GEOS-1, occur in the polar regions. Table 4 shows the global 
mean (rms) regional differences of T s for all ten pairs (ocean and land separately). The largest 
difference appears in (GEOS-1 minus FDTV)forland inJanuary 1992: mean(rms) difference = -4.7 
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(8.6) K. The overall modulus mean difference (i.e., absolute values of the mean differences, not 
affected by subtraction order, see e.g., Reynolds [1988]) (rms difference) is 0.52 (3.3) K (Table 4 
footnotes), very similar to the modulus mean (rms) difference for surface air temperature T a , 0.53 

(2.6) in Zhang et al. [2006]. Like T a , the global bias arises from a partial cancellation between land 
and ocean biases: if land and ocean are separated, the modulus mean (nns) difference of T s is 1.57 

(4.6) and 0.70 (2.4), respectively, compared witb 0.73 (3.7) and 0.77 (1.9) for T a . The land T s 
modulus mean bias is twice that of T a while the ocean values are comparable. 

Figure 8 illustrates the zonal monthly mean differences (all minus FDTV). The largest 
differences, up to 20 K, appear in the polar regions, where cloud contamination makes retrievals 
very difficult. In fact, wintertime polar regimes also exhibit larger clear-sky biases in T s [cf. Curry 
etal., 1996]. However, the reanalyses also exhibit the well-known cold biases at these latitudes. For 
most lower and middle latitudes, the differences are generally ^ 5 K and < 3 K for land and ocean, 
respectively. There is a systematic tendency for the reanalyses to be slightly colder in winter and 
wanner in summer over land relative to FDTV, despite the fact that the original ISCCP T s retrievals 
have clear sky biases (too cold in winter, too warm in summer) [ Rossow and Garden, 1993] that 
have been reduced in the FD analysis [Zhang et al., 2004], 

We have collected three additional SST products: (1) NOAA/NASA Pathfinder Advanced 
Very High Resolution Radiometer (AVHRR) (descending and ascending) version 4.0 with 9 
km/daily resolution and global coverage for 1 98 1 to the current, produced by JPL/Caltech (estimated 
error is given as 0.5-0. 7 K [ Vazquez et al., 1998]); (2) NOAA Optimum Interpolation Sea Surface 
Temperature Analysis version OI.V1 (1983-2000) and V2 (2000-2001) with l°/weekly resolution 
and also global coverage [Reynolds, etal., 2002]; and (3) The Tropical Rainfall Measuring Mission 
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(TRMM) Microwave Image (TMI) version 3 with 25-km/weekly resolution and 40 S - 40 N 
coverage for 1997 to 2004 [Wentz and Meissner, 2000] (also see comparisons in Chelton and Wentz, 
[2005]). The monthly mean time series for these datasets (simplified to “PF” for Pathfinder, “RS” 
for Reynolds, “WZ” for TMI, and “TV” for FDTV as before) are compared (after remapping them 
to the same 280-km equal-area map). 

Table 5 shows global averages of the statistical results from the SST time-series comparison 
for all the individual 280-km cells for each of the six possible pairs of data. The correlation 
coefficients among the first three datasets (Pathfinder, Reynolds and TMI, value rows 1, 2 and 4) 
are all ^ 0.96, despite very different space-time resolutions. The Pathfinder analysis procedure 
actually requires good agreement with the Reynolds dataset. Their mean (rms) difference is ^ 0.3 
(0.3) K, reflecting the fact that all the three datasets are calibrated or regressed against the same in- 
situ SST measurements, primarily buoys, to obtain “bulk” temperatures. The FDTV skin (not bulk) 
temperature values are less well correlated with these three datasets [correlations between 0.69 and 
0.79 with larger mean (rms) differences, ranging from 0.82 (1.1) up to 1.7 (1.7)]. The overall 
modulus mean (rms) difference is 0.71 (0.87), comparable to earlier comparisons in Table 4 (note 
that the reanalyses all use some version of the Reynolds SST). The high bias of the FDTV values 
(by about 1 K) may be accounted for, in part, by the difference between skin and bulk temperatures 
(the former is cooler on average by several tenths of degree at nighttime but can be significantly 
warmer during daytime), but as we show below it is more related to anomalies in the TOVS 
atmospheric temperatures used in the ISCCP retrieval. The larger variations of the FDTV values 
(standard deviation > 2 K compared with ^ 1 .6 K for the other three) is also associated in part with 
the larger variability (especially diurnal) of skin versus bulk temperatures. 
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The zonal-mean differences between these three SST datasets and FDTV for January and 
July, 2000 (not shown) exhibit the same patterns as shown in Figure 8, i.e., the biases between 
FDTV and the reanalyses (and 31) are similar to those with the Reynolds SST, which is used in 
various ways as a standard for all the others. The larger polar differences are all associated with 
increasing sea ice fraction: the other SST products report water temperature, not the ice surface 
temperature as does FDTV. Excluding the polar regions, the difference between FDTV and the other 
datasets is generally ^ 2 K. Figures 9a and 9b also show the differences in the SST anomaly time 
series (all minus FDTV) for global and tropical averages, respectively. Much of the large-amplitude 
variations of the other datasets relative to FDTV are induced in the ISCCP-retrieved T s values by 
spurious temperature/humidity changes in the original TOV S as noted in Zhang et al. [2006: section 
2.1]). 

Because the difference of T s and T a controls the net LW at the surface, Zhang et al. [2006: 
Section 2.9 and Figure 1 1] also compared the differences (T s minus TJ from FDTV with those from 
31 and the three reanalyses but we do not repeat it here. 

3. Summary and Discussion 

Estimating the surface radiation budget has long been pursued [ Simpson , 1929] because, 
together with latent and sensible heat fluxes, it comprises one fundamental coupling of the 
atmosphere with the ocean and land surfaces. Direct estimates that resolve regional and weather- 
scale variability with reasonable accuracy have only become possible with the advent of complete 
global, mostly satellite, datasets within the past couple of decades. Zhang et al. [2004] now estimate 
that surface radiative fluxes can be determined to within about 10-15 W/nr (regional monthly mean 
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flux values), improved by about 5 W/m 2 over Rossow and Zhang \ 1995]. The main limitation is still 
the accuracy of the input datasets. Sensitivity studies using ISCCP -based-estimates (“realistically- 
assumed”) of the input uncertainties [based, in part, on the analysis presented here and in the work 
of Zhang et al., 1995; 2004; 2006] show that the leading uncertainties in the surface fluxes are no 
longer predominantly induced by clouds but are now as much associated with uncertainties in the 
surface and near-surface atmospheric properties. 

This study presents a fuller, more quantitative evaluation of these uncertainties for the 
surface albedo and emissivity, and surface skin temperatures by comparing the main available global 
datasets that are treated as an ensemble of realizations of the actual climate such that their 
differences represent an estimate of the uncertainty in their values because we do not know the 
“truth”. The results are globally representative and may be taken as a generalization of our previous 
ISCCP-based estimates in Zhang et al. [1995; 2004]. Generally speaking, the radiative flux 
uncertainties induced by the input uncertainties from Zhang et al. [2006] and this study are 
consistent with the flux uncertainty estimates in Zhang et al. [2004] based on all of our past and 
current sensitivity studies, and therefore, in order to meet the accuracy requirements for 
climatological studies, say within 10 W/nr, we still need substantial improvements of those input 
parameters that contribute maj or errors as summarized in the work of Zhang et al. [2006] and below. 

The main conclusions ofthis study, together with all the radiative flux sensitivity studies we 
have done previously, are as follows. (1) The uncertainty in global mean (rms) for regional 
differences for visible, NIR and total SW black-sky albedo (land only) is 0.3% (7.5%), 16% (9%) 
and 7.6% (7%), respectively. Although land surface albedos in the NIR remainpoorly constrained, 
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they do not cause too much error in SW fluxes because the atmosphere absorbs much of the 
radiation in this wavelength range. Nevertheless, the subtle geographic and seasonal variations of 
surface absorbed SW depend on variations of the NIR albedo influenced by the precise mixture of 
soil/rock, vegetation and snow and the seasonal variations of vegetation activity and snow. The 
mean (nns) difference of 7.6% (7%) in broadband black-sky SW albedo disagreement for land 
surface from this study leads to the uncertainty estimate of about 7%, which can easily induce 5-10 
W/m 2 uncertainty in (up welling) surface SW flux estimates. But the full potential of satellite sensors 
like MODIS or the MEdium Resolution Imaging Spectrometer Instrument (MERIS) that sample 
(nearly) the complete solar wavelength range to produce definitive surface albedo maps with 
relatively high space-time resolution has not yet been achieved. The current 16-day period MODIS 
albedo, though having very high spatial resolution, is usually not globally complete. The current 
product is also produced with a 16-day time resolution. Most importantly, the current product 
corrects for atmospheric effects with a relatively crude climatology of atmospheric composition and 
aerosols: neglecting the detailed space-time variations of water vapor, in particular, makes the 
uncertainties for the NIR albedo large. These limitations are not inherent; given the available data 
products specifying the ozone, water vapor and aerosol abundances as functions of location and 
time, a much more intensive analysis of the MODIS radiances is possible. However, to obtain the 
surface radiative fluxes over the whole time period of the ISCCP-FD product, the MODIS product 
must be used to “calibrate” the longer AVHRR record for surface albedo. (2) In the same way, the 
uncertainty for broadband and window emissivity may be taken as about 3-5%. Even though the 
disagreement between the two surface LW emissivity reconstructions appear to be relatively large, 
the fact that this disagreement is confined to wavelengths > 20 pm means that it has little practical 
effect (1-3 W/m 2 ) on the surface upwelling LW fluxes because the atmosphere is nearly opaque at 
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these wavelengths, so reflected downwelling radiation almost balances the reduction of upwelling 
radiation (this near-balance is not true in the polar regions but there the fluxes are much smaller in 
magnitude). (3) The uncertainty in the skin temperatures is about 3 K. It is one of the two leading 
factors that cause problems with surface LW fluxes. Even though the disagreements among the 
various datasets are generally only 2-4 K, this can easily cause 10-15 W/nr uncertainty in surface 
(upwelling) LW calculation. 

The overall summary, based on both Zhang et al. [2006] and this work, is as follows. 

For surface LW fluxes, significant improvements could be obtained by improving the retrievals of 
(in order of decreasing importance): (1) surface skin temperature, (2) surface air and near-surface- 
layer temperature, (3) column precipitable water amount and (4) broadband emissivity. As we have 
shown the retrieval of surface skin and air temperatures are linked because the relevant satellite 
measurements are sensitive to both so they should be retrieved together [cf. Aires et al., 2002]. 
Moreover, the systematic differences between clear- and cloudy-sky skin temperatures requires 
extension of the analysis of infrared measurements to include microwave measurements [cf. Prigent 
etal., 2003]. Some improvement in surface emissivities over a limited wavelength range may come 
from a careful analysis of the new infrared spectrometers, such as the Atmospheric Infrared Sounder 
[AIRS, see Olsen, 2005], For surface SW fluxes, improvements could be obtained (excluding 
improved cloud treatment) by inproving the retrievals of (1) aerosols (from our sensitivity studies 
but not discussed in this work), and (2) surface (black-sky) albedo, of which, NIR has larger 
uncertainty. Surface albedo and aerosols are also linked since satellite measurements at solar 
wavelengths are generally sensitive to both, especially over land, suggesting again, that the best 
results would be obtained by retrieving both together. 
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FIGURE CAPTIONS 


Figure la. Scatter plot of visible black-sky albedo (in %) for Amazon Rainforests (dark grey 
diamond), snow-free Grassland (empty circle) and snow-free deciduous forests (solid triangle) 
from ISCCP6 and MODIS for four 16-day periods (beginning on April 6, July 1 1 and September 
29, 2000, and January 17, 2001). 

Figure lb. Same as la but for near-infrared and the diagonal line is not plotted from origin (0,0), 
but from the point that visually approximately reflect their mean bias for all the corresponding 
regions. 


Figure lc. Same as lb but for total SW (visible and near-infrared). 

Figure Id. Scatter plot of visible black-sky albedo (in %) for snow-covered grassland (empty 
triangle), land ice (dark grey square), sea ice (solid triangle) and Sahara deserts (dark grey 
diamond) from ISCCP6 and MODIS for four 16-day periods (beginning on April 6, July 1 1 and 
September 29, 2000, and January 17, 2001). 

Figure le. Same as Id but for near-infrared and the diagonal line is not plotted from origin (0,0), 
but from the point that visually approximately reflect their mean bias for all the corresponding 
regions. 

Figure If. Same as le but for total SW (visible and near- infrared). 

Figure 2a. Global map of differences (ISCCP6 minus MODIS) of visible black-sky albedo (in 
%) from four 16-day periods (beginning on April 6, July 1 1 and September 29, 2000, and 
January 17, 2001). 

Figure 2b. Same as 2a but for near-infrared. 

Figure 2c. Same as 2a but for total SW (visible and near-infrared). 

Figure 3a. Scatter plot of monthly-mean TOA clear-sky SW-up fluxes (in W/m 2 ) for snow-free 
land from ISCCP-FD and CERES, accumulated from January and July 1992. 

Figure 3b. Same as figure 3a but for MODIS-FD and CERES. 

Figure 3c. Same as figure 3a but for snow-covered land. 

Figure 3d. Same as figure 3b but for snow-covered land. 

Figure 4a. Comparison of zonal monthly-mean TOA clear-sky albedo (in %) from ISCCP-FD, 
SRB and ERBE for July 1986. 
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Figure 4b. Same as figure 4a but for January 1987. 

Figure 5. Global map of differences (ISCCP-FD minus SRB) of broadband emissivities (in %) 
from ISCCP-FD 1992 annual mean and SRB climatology. 

Figure 6a. Water emissivity from theoretical calculation using NASA GISS model (see text) at 
wind speed = 3 m/s for cosine zenith view angle = 0.1, 0.5 and 1.0, respectively. 

Figure 6b. Same as figure 6a but for a fixed cosine zenith view angle = 0.5 and for wind speed = 
3 and 7 m/s, respectively. 

Figure 7. Global map of differences (ISCCP-FD minus SRB) of window-region emissivities (in 
%) from ISCCP-FD 1992 annual mean and SRB climatology. 

Figure 8a. Difference of zonal monthly-mean surface skin temperature (in K) for ERA15, 
GEOS-1, NCEP and 31 minus FDTV, respectively (as TSOlertv for ER - TV, etc.), for land and 
January 1992. 

Figure 8b. Same as figure 8a but for July 1992. 

Figure 8c. Same as figure 8a but for ocean. 

Figure 8d. Same as figure 8b but for ocean 

Figure 9a. Difference of the time series of global mean SST (in K) from Pathfinder (PF), 
Reynolds (RS) and TMI (WZ) minus FDTV (TV), respectively. 

Figure 9b. Same as figure 9a but for 30 S - 30 N mean. 
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TABLES 


Table la. Comparison between ISCCP6 and SRB for monthly-mean broadband clear-sky surface 
albedo 3 for snow-free land, ice-free water, snow-covered land and sea ice b for January, April, July, 
October, 1992: globally-averaged mean difference (ISCCP6 minus SRB) with its spatial standard 
deviation in parentheses from 280-km equal-area maps. Last row shows the mean values of the 
above four months. 



Snow-free Land 

Ice-free Water 

Snow-covered Land 

Sea Ice 

9201 

-1.7 

( 5.8) 

-0.7 

(1.9) 

15.0 i 

! 11.4) 

20.6 

( 16.8) 

9204 

-1 . 0 

( 4.5) 

-0 . 3 

( 3.8) 

10.9 ' 

( 3.5) 

7 . 5 

( 6.4) 

9207 

-3 . 4 

( 3.8) 

-0 . 4 

( 4.4) 

7.1 

( 4.8) 

8 . 1 

( 3.8) 

9210 

-2 . 0 

( 4.4) 

0 . 0 

( 2.1) 

11.7 

( 6.3) 

7 . 8 

( 6.1) 

mean 

-2 . 0 

(4.6) 

-0.3 

( 3.0) 

11.2 

( 6.5) 

11.0 

( 8.3) 


“In % 


b Defined as > 95% of a grid box through all the days of a month using ISCCP snow and ice data. 


Table lb. Comparison between ISCCP6 and SRB for monthly-mean broadband clear-sky surface 
albedo a , collected from the four months (Table la) for seven regions 11 


region 

number 

ISCCP6 
(X) mean 

SRB 

(Y) mean 

X-Y 

mean 

Stdv of 
(X-Y) 

Regression 

grid box 
number 

corr. Coef 

slope 

intercept 

Nm Dev 

i 

14.5 

19.9 

-5 . 4 

1 . 9 

0.714 

1.4 

0 . 1 

0 . 8 

160 

2 

32 . 8 

29.9 

2 . 9 

4 . 3 

0.766 

0.8 

3 . 6 

3 . 0 

484 

3 

15.2 

16.2 

-1.1 

3 . 5 

0.515 

0.8 

3 . 5 

2 . 4 

150 

4 

11.6 

14.0 

-2 . 4 

3 . 1 

0.607 

0.9 

3 . 6 

2 . 2 

277 

5 

54.5 

44.0 

10.5 

4 . 9 

0.973 

0.8 

4 . 8 

2 . 5 

33 

6 

76.6 ' 

73.2 ’ 

3 . 4 

2 . 6 

-0.280 

0.9 

2 . 2 

2 . 1 

152 

i 

6 9.0 

61.0 

8 . 0 

6 . 1 

0.773 

0.8 

4 . 8 

3 . 5 

T74 


a In % 


b As numbered in column one, the seven types of land features (based on ISCCP surface information) are: (1 ) Amazon 
rainforests (2) Sahara Deserts, (3) Snow-free grassland, (4) Snow-free deciduous forests, (5) Snow-covered grassland, 
(6) Land ice (w ith or without snow cover), and (7) S ea ice. All surface features (other than vegetation types) are defined 
when such features > 95% of a 280-km equal-area box and all vegetation types are defined as they are > 66% of a grid 
box. “Nm Dev” is the RM S distance o f all points from the regression line, “corr.co ef ’ correlation coefficient, and “stdv” 
standard deviation. 
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Table lc. Comparison between ISCCP6 and MODIS for 16-day mean visible, near-infrared (NIR) 
and total SW broadband black-sky land surface albedo (in %) from four 16-day periods (see text) 
for: (1) snow-free land, (2) snow-covered land 3 , and (3) all land on 280-km equal-area maps. 


Band 

ISCCP6 

MODIS 

X-Y 

Stdv of 

Regression 

grid box 


(X) mean 

(Y) mean 

mean 

(X-Y) 

corr. Coef 

slope 

intercept 

Nm Dev 

number 



(1) Snow-free Land (when both are 

snow-free) 



visible 

11.9 

11.0 

0.87 

3.53 

0.9093 

0.995 

-0.81 

2.50 

3842 

NIR 

11.8 

30.9 

-19.12 

6.69 

0.7831 

0.998 

19.14 

4.74 

3842 

Ttl SW 

11.9 

21.0 

-9.13 

4.32 

0.8771 

0.993 

9.21 

3.07 

3842 

(2) 

Snow-covered Land (when both 

snow >75%) 



visible 

44.3 

47.3 

-2.99 

14.25 

0.8584 

0.995 

3.22 

10.11 

1137 

NIR 

31.4 

38.2 

-6.81 

10.09 

0.7626 

0.742 

14.91 

7.50 

1137 

Ttl SW 

38.5 

42.7 

-4.18 

11.63 

0.8400 

0.902 

7.95 

8.52 

1137 




(3) All Land (when both data are ; 

available) 




visible 

18.3 

18.0 

0.31 

7.51 

0.9318 

1.058 

-1.37 

5.11 

5958 

NIR 

15.6 

31.4 

-15.74 

9.05 

0.7259 

0.683 

20.69 

6.71 

5958 

Ttl SW 

17.2 

24.7 

-7.55 

6.96 

0.8957 

0.880 

9.61 

5.04 

5958 


“Defined as both are snow-free and >75% snow of a ce' 


1, respectively. 


Table Id. Comparison of monthly mean surface clear-sky albedo (%) and upwelling SW flux (Wm 
2 ) beteen ISCCP-FD andMODIS-FD (see text) for collective July 2000 and January 2001 monthly 
means for: (1) snow-free land, (2) snow-covered land, and (3) all land on280-km equal-area maps. 


Clear-sky 

FD 

(X) mean 

MODIS-FD 
(Y) mean 

X-Y 

mean 

Stdv of 
(X-Y) 

Regression 

grid box 
number 

corr. Coef slope 

intercept Nm Dev 

(1) Snow-free Land (when both are snow-free) 

Albedo 

14.4 

15.7 

-1.33 

3.50 

0.9287 

0.868 

3.23 

2.47 

1981 

SW-up 

40.4 

44.3 

-3.87 

10.13 

0.9218 

0.934 

6.56 

7.30 

1981 

(2) Snow-covered Land (when both snow >75%) 

Albedo 

49.5 

45.2 

4.32 

13.26 

0.8166 

0.800 

5.56 

9.76 

568 

SW-up 

74.1 

74.5 

-0.69 

12.51 

0.9941 

1.052 

-3.15 

7.85 

568 

(3) All Land (when both data are available) 

Albedo 

21.1 

21.0 

0.14 

7.12 

0.9249 

0.832 

3.41 

4.91 

3002 

SW-up 

46.4 

48.9 

-2.51 

11.09 

0.9775 

1.017 

1.70 

7.75 

3002 
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Table 2a. Same as Table lb but comparing ISCCP-FD and ERBE for broadband TOA clear-sky 


albedo (in %) collected from monthly means for January, April, July and October of 1986. 


region 

ISCCP-FD 

ERBE 

X-Y 

Stdv of 

Regression 

grid box 

number 

(X) mean 

(Y) mean 

mean 

(X-Y) 

corr. Coef 

slope 

intercept 

Nm Dev 

number 

i 

14.5 

15.6 

-i.i 

1 . 4 

0.386 

1 . 0 

1 . 5 


124 

2 

31.3 

32 . 3 

-1.0 

2 . 2 

0.866 

1 . 0 

1 . 3 

1 . 5 

484 

3 

18.4 

19.2 

-0.8 

3 . 0 

0.450 

1 . 0 

1 . 8 

1 . 7 

150 

4 

18.0 

17.6 

0 . 4 

3 . 2 

0.241 

1 . 0 

-0.1 

2 . 0 

272 

5 

43.0 

40.1 

2 . 9 

12.9 

0.740 

1 . 0 

0 . 4 

3 . 3 

40 

6 

65 . 7 

73.5 

-7.8 

2 . 5 

0.032 

1 . 1 

-2 . 4 

1 . 8 

148 

7 

57 . 8 

63 . 3 

-5.5 

6 . 3 

0.447 

1 . 1 

-1.5 

2 . 8 

133 


Table 2b. Same as Table 2a but for SRB and ERBE. 


region 

SRB 

ERBE 

X-Y 

Stdv of 

Regression 

grid box 

number 

(X) mean 

(Y) mean 

mean 

(X-Y) 

corr. Coef 

slope 

intercept 

Nm Dev 

number 

i 

17.6 

15.6 

2 . 0 

1 . 3 

0.339 

1 . 2 

-4.2 

2 . 0 

124 

2 

30.0 

32 . 3 

-2 . 2 

2 . 5 

0.821 

1 . 0 

2 . 0 

1 . 8 

484 

3 

19.3 

19.2 

0 . 0 

2 . 6 

0.309 

1 . 0 

0 . 7 

1 . 7 

150 

4 

18.0 

17.6 

0 . 4 

2 . 5 

0.129 

1 . 0 

-0.5 

1 . 7 

272 

5 

36.5 

40.1 

-3 . 7 

13.0 

0.766 

1 . 1 

-2 . 1 

2 . 9 

40 

6 

60.2 

73.4 

-13.2 

2 . 1 

0.512 

1 . 3 

-5 . 2 

1 . 9 

144 

7 

52 . 5 

61.4 

-8.9 

7 . 0 

0.204 

1 . 2 

-3 . 1 

2 . 8 

110 


Table 2c. Comparison between ISCCP-FD and CERES for monthly-mean TOA clear-sky albedo 
collected from July 2000 and January 2001 (snow maps are from FD and MODIS only). 


ISCCP-FD 

CERES 

X-Y 

Stdv of 

Regression 

grid box 

(X) mean 

( Y ) mean 

mean 

(X-Y) 

corr. Coef 

slope 

intercept 

Nm Dev 

number 

(1) Snow-free Land 

| is . ! ; 

18.4 

| - 0 . 2 9 

| 2 . 38 | 

| 0 . 9231 | 

1 0 . 9 1 3 | 

! . 87 

! . 7! j 

1 1 9 7 8 ] 

I < 2 

) Snow-covered Land [ 

| 44 . 7 ' 

| 40 . 0 

| 4 . 64 

] 11.34 | 

| 0 . 8082 | 

1 1 . 0 8 4 | 

| -8.42 

| 7 . 65 | 

1 48 6 i 

j (3) All Land j 

| 22 . 8 ; 

| 2! . 9 

| 0 . 90 

| 5 . 70 ; 

|0.8972 j 

1 0 . 8 8 0 | 

1 1.84 

| 4 . 12 j 

I 2916 1 


Table 2d. Same as Table 2c but for MODIS-FD and CERES. 


MODIS-FD 

CERES 

X-Y 

Stdv of 

Regression 

grid box 

(X) mean 

(Y) mean 

mean 

(X-Y) 

corr. Coef 

slope 

intercept 

Nm Dev 

number 

(1) Snow-free Land 

| 19 . 1 

18.4 

I 0 . 70 


| 0.9616 | 

mm,: 

| -0.87 

| 1 . 16 j 

1978 ( 

< 2 

) Snow-covered Land | 

| 43 . 0 ; 

| 40 . 0 

2.97 

9.76 

| 0.8703 j 

1.175 

| - 1 0 . 5 2 

| 6 . H j 

| 486 I 

j (3) All Land | 

| 23.1 

| 21.9 | 

| ! . 21 

| 4 . 50 

| 0.9323 | 


| -0.94 

| 3.20 | 

[2916 1 
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Table 3a. Comparison of monthly-mean clear-sky broadband planetary albedo (in %) and clear-sky 
TOA upwelling SW flux (in W/nr) between the results from four different radiative transfer models 
with their own surface albedos (ISCCP2-ALB, GCM-ALB, ISCCP6-ALB and SRB-ALB) as well 
as other inputs and the ERBE results for July 1986 and January 1987 a 



a The statistical numbers are globally-averaged mean differences (all minus ERBE) with their corresponding spatial 
standard deviations in parentheses. The radiation models used to produce the results are old NASA GISS GCM [as 
ISCCP2, see Zhang et al., 1995], current GISS model (as GCM), revised current GISS model [as ISCCP6 = 03-model 
in Zhang et al., 2004], and WCRP/GEWE X SRB model (release 2) [as SRB, see Stackhouse et al., 1999], respectively. 
Their main input datasets are ISCCP-C 1 forISCCP2, andlSCCP-Dl for the rest three models. Snow-free land and ice- 
free ocean areas are separated, defined as such surface features are > 95% of a 2.5° equal-area box through all the days 
of a month using ISCCP6. 


Table 3b. Same as Table 3a, except for snow-covered land and sea ice, respectively. 






















































































Table 4. Comparison of monthly mean surface skin temperature from ERA 15, GEOS-1, NCEP, 
FDTV and 31 for the four seasonal months of 1992 a 


X/Y 

Jan. 1992 

April 1992 

July 1992 

Oct. 1992 

Ocean 

Land 

Ocean 

Land 

Ocean 

Land 

Ocean 

Land 

ER/GE 

0.57 ( 2 . 8 ) 

1.21 ( 4 . 9 ) 

0.33 ( 2 . 1 ) 

- 0.23 ( 3 . 8 ) 

0.67 ( 3 . 4 ) 

- 0.27 ( 3 . 3 ) 

0.70 ( 3 . 0 ) 

0.94 ( 4 . 2 ) 

ER/NC 



0.12 ( 1 . 2 ) 


0.24 ( 1 . 7 ) 

- 0.27 ( 3 . 5 ) 

0.33 ( 1 . 5 ) 

- 0.66 ( 3 . 6 ) 

ER/TV 

1.22 ( 2 . 7 ) 


1.68 ( 2 . 3 ) 

- 3.25 ( 5 . 3 ) 

2.05 ( 2 . 1 ) 

0.23 ( 5 . 4 ) 

1.12 ( 2 . 7 ) 

- 2.89 ( 4 . 7 ) 

ER/3I 

- 0.11 ( 1 . 8 ) 


0.27 ( 1 . 2 ) 


0.67 ( 2 . 2 ) 

- 0.79 ( 4 . 6 ) 

0.37 ( 2 . 3 ) 

- 1.06 ( 4 . 0 ) 

GE/NC 


- 3.93 ( 5 . 9 ) 

- 0.21 ( 1 . 6 ) 


- 0.43 ( 2 . 3 ) 

0.00 ( 3 . 0 ) 

- 0.37 ( 2 . 1 ) 

- 1.61 ( 3 . 9 ) 

GE/TV 

0.64 ( 4 . 6 ) 


1.36 ( 3 . 1 ) 


1.38 ( 4 . 0 ) 

0.50 ( 5 . 2 ) 

0.43 ( 3 . 7 ) 

- 3.83 ( 5 . 0 ) 

GE/3I 

- 0.68 ( 3 . 3 ) 


- 0.06 ( 2 . 0 ) 


0.00 ( 2 . 3 ) 

- 0.84 ( 4 . 5 ) 

- 0.33 ( 2 . 7 ) 

- 2.33 ( 4 . 7 ) 

NC/TV 

1.24 ( 2 . 6 ) 

- 0.81 ( 5 . 4 ) 

1.56 ( 2 . 4 ) 

- 1.77 ( 4 . 2 ) 

1.81 ( 2 . 6 ) 

0.49 ( 4 . 3 ) 

0.80 ( 2 . 5 ) 

- 2.22 ( 4 . 6 ) 

NC/3I 

- 0.08 ( 1 . 8 ) 

0.85 ( 3 . 6 ) 

0.15 ( 1 . 3 ) 

- 0.09 ( 3 . 9 ) 

0.43 ( 1 . 6 ) 

- 0.95 ( 4 . 0 ) 

0.04 ( 1 . 5 ) 

- 0.87 ( 3 . 3 ) 

TV/31 

- 1.32 ( 2 . 5 ) 

1.95 ( 5 . 5 ) 


1.35 ( 4 . 8 ) 



- 0.75 ( 2 . 5 ) 

1.48 ( 4 . 8 ) 

abs mean 

0.65 ( 2 . 6 ) 

2.42 ( 5 . 6 ) 

0.72 ( 1 . 9 ) 

1.50 ( 4 . 3 ) 

0.91 ( 2 . 5 ) 

0.59 ( 4 . 4 ) 

0.52 ( 2 . 5 ) 

1.79 ( 4 . 3 ) 


tl The shown numbers are global land/oceanic means and standard deviations of regional differences (X - Y) on 280-km 
equal-area grid map in Kelvin. Ocean areas include coast cells (defined as area> 33% of water; otherwise, land). “ER”, 
“GE”, “NC”, “TV” and “3P areforERA15, GEOS-1 , NCEP , FDTV and 31, respectively. The ten X/Y pairs are all their 
possible combinations. The last row (“abs mean”) is the modulus (i.e., absolute-value) average ofthe mean differences 
(standard deviation) ofthe ten pairs, and it is not affected by the X/Y assignment. When averaging over all the four 
months of the last row, the pseudo-annual/all-pair mean-moduli (rms) differences = 0.70 (2.4) and 1 .57(4.6) K, for ocean 
and land, respectively. F or the global (ocean and land) mean, it becomes 0.52 (3.3 ) K. 
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Table 5. Global mean statistical results from the SST time series comparison for all the individual 
280-km grid cell for all the possible six pairs from four datasets 3 


X/Y 

X 

Y 

M Diff 

Stdv 

Cor. Coef 

Slope 

Intercept 

Nm Dev 

PF/RS 

292.2 (1.7) 

292.4 (1.6) 

-0.28 

0.32 

0.958 

0.913 

25.54 

0.22 

PF/WZ 

297.5 (1.6) 

297.8 (1.6) 

-0.28 

0.29 

0.973 

0.992 

2.79 

0.20 

PF/TV 

292.2 (1.7) 

291.3 (2.6) 

0.82 

1.60 

0.712 

1.168 

-48.13 

1.00 

RSAVZ 

297.7 (1.6) 

297.8 (1.6) 

-0.07 

0.25 

0.981 

1.033 

-9.93 

0.16 

RS/TV 

292.2 (1.6) 

291.1 (2.7) 

1.11 

1.70 

0.689 

1.217 

-61.00 

1.02 

WZ/TV 

297.8 (1.6) 

296.0 (2.0) 

1.72 

1.06 

0.786 

1.047 

-15.75 

0.70 


“As described in the text, the four datasets, Pathfinder, Renolds, TMI and FDTV, are denoted asPF, RS, WZ and TV 
in the table. The statistical values are based on the coordinates of points in a scatterplot with Y values on the ordinate 
and X values on the abscissa. Regression statistics are from a linear least squares fit to the scatter of points. All values 
are in K, except the correlation coefficients and slopes, which are unitless. The first two value columns are mean (spatial 
standard deviation) for X and Y, respectively. “M Diff” (“stdv”)is mean (rms) difference between X and Y, and "Nm 
Dev” the rms distance of all the points from the regression line. The statistical results are from the compared time series 
for all the individual 280-km cells over the globe. All the datasets have a global coverage except WZ that only covers 
40 S to 40 N ocean. The matched time periods vary from pair to pair with a maximum coverage from July 1983 to 
December 2003 (see text). 
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Fig. 3a 
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TOA clear-sky SW-up from ISCCP-FD (W/m A 2) 
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Fig. 3b 
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TOA clear-sky SW-up from MODIS (W/m A 2) 
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Fig. 3c 
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Fig. 3d 
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Fig. 9a 
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Fig. 9b 
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